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INPUT DESCRIPTION FOR JAMESON'S THREE-DIMENSIONAL 
TRANSONIC AIRFOIL ANALYSIS PROGRAM 
By Perry A. Newman aad RnLy-M.,- Davis 
Langley Research Center 


SUMMARY 


This paper describes the input parameters for a computer program which 
has been developed by Antony Jameson of the Courant Institute of Mathematical 
Sciences, New York University under NASA Grant NGR-33-0l6-l67 . References 
1 and 2 describe the method and give several sample restilts; however, a user’s 
manual (ref, 3), is not yet available. A version of the program which was 
successfully demonstrated on the CDC system here at NASA Langley is available 
as Langley Program A^f231. This paper gives a description of the input 
parameters and a listing of input data cards and some output results for 
a sample case. It should enable one to use this program until a final version 
with its formal user’s manual is available. 

A brief description of the program and type of problems it will handle is 
given in the computer program abstract which is reproduced as Figure A1 of 
Appendix A. Further details concerning the method are given in references 1 
and 2. An overlay diagram and a list of subroutine names and functions are 
given as Figures P2 and A3 respectively. 

Computer storage and time requirements for a three-dimensional problem 
limit what can reasonably be done on present computers. Therefore, in 
practice, this 3-D program does not have the flexibility (with respect to 
grid size, number of grid halvings, and multiple case runs) of similar 
2-D analysis programs (refs. U-7) which are based on the full nonlinear 
potential formulation. In this program, the potential function G is a 
large three-dimensional array which, in essence, determines both the central 
memory storage and the computational time. Typical runs, at the dimensions 
stored on the data cell, require (a) 15-30 minutes CPU time on a CDC 6600 
(using a run compiler) to reach a convergence criterion of about 10“5 and 
(b) 322K (base 8) central memory storage. These runs would start on a 
coarse grid {6h x 8 x l6) and be halved only once. Iterations take very 
close to 8 times as long on the resiilting fine grid (128 x l6 x 32). 

The maximum values quoted here, and on the following pages are those consistent 
with the data cell dimensioning and one halving of the grid; these , of 
course, can be changed somewhat. In practice there is little advantage in 
using the multiple case capability; thus it is recommended that single cases 
be run. 
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INPUT DESCRIPTION 


Description and Comments 


TITLE. - Descriptive title of case or sequence; 
Format (8 a10) 

Appears on Varian Plots and "beginning 
of output, 

DESC. - Description for card in Read Order 3 
Format (8 a10). 

FNX, PNY, FNZ, FPLOT 

Format (8E10.7) 

Note; A number of quantities are read in 
as floating-point numbers and converted 
to integers vithin the program. 

FNX. - Number of computational grid points in 

"chordvise direction” from downstream 
infinity* soround the leading edge and 
back to downstream infinity on coarsest 
mesh. Maximum is 64 (128 with no grid 
halving). 

FNY. - Number of coaqputational grid points in 

'’normal direction" from airfoil surface 
to infinity on coarsest mesh. Maximum is 
8 (l6 with no grid halving). 

FNZ. - Number of computational grid points in 

"spanwise direction" from infinity, 
across the wing span and to infinity on 
coarsest mesh. Maximum is l6 (52 with no 
I grid halving ) . 


PPLOT. - Plot trigger. Selects type of plot for 
chordwise surface pressiire coefficients. 
PPLOT = 0. Printer plots, one at each 
spanwise grid plane section with CP 
versus the conqputational grid chordwise 
variable. 



Head 

Order 


Number 

Cards 


Description and Comments 

FPLOT =! 1* Varlan plots (from !EHREED)» 
These are superimposed plots^ vlth all 
span sections shown on two figures, an 
upper surface and a lower surface plot 
of CP versus physical space chord- 
wise variable. 

FPLOT =2. Varlan plots (from THREED) 
as above plus section plots (from GRAPH), 
These latter plots, one per section, 
give upper and lower siirface CP versus 
physical space chordwise variable. 


Defaults to zero 


1 

DESC. 

Description for card in Read Order 5 
Format (8A10). 

1 card 
for each 
computational 
grid. Maximum 
essentially 2, 
dimensioned 5. 

FIT, COVO, 

FIT. 

COVO. 

PIO, P20, P3Q, BETAO, STRIPO, FHALF 
Format (8E10',7) 

Maximum number of iterations on 
this gird, called ffiT in program. 

Convergence criterion on the mayiiTinTn 
change in reduced velocity potential 
(G) from one iteration cycle' to the 
next on this grid. , 

1 


PIO. 

Subsonic point relaxation factor on 
this grid; must be <2.. Typically 
1.6 on coarse grid. 


P20. 

Supersonic point relaxation factor; 
must be < 1.. Should use 1.0 for- 



stability on all grids. 


P30. 

Circulation relaxation factor. May 
be > 1.0 . 


BETAO. -■ 

Stabilization factor used at super- 
sonic points in finite difference 
operator if BETAO > 0, Most needed 
when M^ >1. , many cases operate satis- 
factorily with BETAO = 0. Convergence 
is slowed but stability enhanced when 
BETAO >0. 
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Read 

Order 


Number 

Cards 


Description and Comments 


6 1 

7 1 


STRIPO. - 

Line relaxation control. Computational 
X-Y planes are relaxed by horizontal 
lines (YSWEEP) in central strip, verti- 
cal lines (XSWEEP) in outer strips. 
STRIPO specifies the fraction of compu- 
tational plane included in central 
strip: 0. _< STRIPO < l1 , where STRIP0=1. 

gives all horizontal line relaxation. 

FHALF. - 

Grid halving trigger. 

1 FHALF 1^1. read another card (Read 
Order 5 format) containing computation- 
al parameters to be used on grid with 
mesh size halved in all directions. 
|fHALF|< 1. must appear on finest grid 
card (last one read). Calculation pro- 
ceeds automatically through the sequence 
of computational grids. 

DESC, 

Description for card in Read Order 7 
Format (8 a10). 

FMACH, YA, 

AL, CDO 

Format (8E10.7) 

FMACH. - 

Freestream Mach number. 

YA. 

Yaw angle (in degrees). 

AL. 

Angle of attack (in degrees) measured 
in plane normal to leading edge, not in 
plane containing freestream direction. 

CDO. 

Drag coefficient due to skin friction 
(CD FRICTION on output). This input 
number is added to the drag coefficient 
obtained by integrating the surface 
pressures (CD FORM on output). 


Read Orders 8 throu^ 19 are used to specif the wing geometry (in physical 
space, of course). One can define the wing at up to 11 span stations. A 
set of airfoil coordinates must he read in at the first station. It need not 
he read in at other stations , if one is changing only combinations of the 
following three airfoil section parameters: chord, thickness ratio or 

^gle of attack (twist). The wing shape at intermediate span positions 
(i.e., the computational grid planes for example) is obtained hy linear 
interpolation in the spanwise direction in the physical space. 
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Read 

Order 


Number 

Ceirds 


Description and Comments 


A multiple run capability for the .same wing geometry at several flow conditions 
is available and controlled by the parameter FWC in Read Order 9. One does 
not have to re-read the wing geometry to use it. The following description 
of Read Orders 8 through 19 is that: required for a single case. 


Read Orders 8 ^d 9 are read only once; 10 and 11 are read FNC (see 9) times; 
12 though 17 (19 if non-symmetric airfoil section) must be read at first 
section and may be required at other sections, depending on the wing geometry. 


8 

1 

DESC. 

- Description for card in Read Order 9 
Format (8A10). 

( 

9 

1 

5NC 

i 

Format (8E10.7) 



FNC. 

Number of span stations at which the 
wing is described or specified. Maximum 
is 11. Must be at least 2 for a 
single case or the first of a sequence. 
If PNC < 2. geometric wing data is 
assumed to be the same as in previous 
case and calculation begins for the new 
flow conditions reading _no further 
input cards. 

10, 

1 

DESC. 

- Description for cards in Read Order 11 
Format (8A10). 

11 

1 

ZS(K), 

CHORD, THICK. AL, FSEC 
Format (QeIO^'J) 



ZS(K). 

— Spanwise coordinate of the wing section 
being specified. It is in the same 
units as CHORD. These stations are 
ordered from tip-to-tip, in ascending 
algebraic order of ZS(K). 



CHORD, 

- Section chord length. The chord of 
the airfoil coordinates to be, read in 
(or already read in at the prior 
station) will he scaled to this 
value . 



■ THICK. 

- Section thickness ratio relative to 
that of the airfoil coordinates to be 
read in (or already read in at the 
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Read 

Order 


Number 

Cards 


Description and CoTnment 


prior station). Note, this is a ratio 
of thichness/chord ratios. The thick- 
ness of the airfoil coordinates will 
be scaled with this value, 

AD. - Section angle of attack or twist 

(in degrees). Airfoil coordinates will 
be rotated through this angle. 

FSEC. - Section airfoil coordinate trigger. 

FSEC =0. Do not read airfoil coor- 
dinates. Last set of airfoil coordin- 
ates read will be used at this section. 
They may be scaled by any combination 
of CHORD, THICK, or AL read above. 
Skip Read Orders 12 through 19 for 
this section. 

FSEC = 1. Read a new set of airfoil 
coordinates which will be used at this 
station and perhaps at other stations. 
They may be scaled by any combination 
of CHORD, THICK, or AL read above 
for this section. 

At first station (K = 1) FSEC is 
ignored; one must supply Read Orders 
12 through 17. 

12 1 DESC. - Description for cards in Read Order 13 

Format (8 a10). 

13 1 FSYM, FNU, FNL 

Format (8E10.7) 

FSYM. - Airfoil symmetry trigger. 

FSYM ^ 1, Symmetric sdrfoil. 

Read in only upper surface airfoil 
coordinates, ordered leading edge 
to trailing edge. 

FSYM <1. Non-symmetric airfoil. 

Read in upper and lower stirface 
airfoil coordinates, respectively, 
each set ordered leading edge to 
trailing edge. Note that leading- 
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Read 

Order 

Number 

Cards 


Description and Comments 




edge points are 
included in both sets. 



PNU. 

Number of coordinates read in for 
upper surface of airfoil. 



FNL. 

Number of coordinates read in for 
lower surface of airfoil. 

Ik 

1 

DESC. 

Description for cards in Read- Order 15 
Format (8 a10). 

15 

1 

TRL, SLT, 

XSING, YSING 

Format ( 8 EIO. 7 ) ' 



TRL. 

Included angle of trailing edge of 
airfoil (in degrees). 

1 



SLT. 

Slope of airfoil mean camber line at 
trailing edge. 



XSING. - 

X coordinate of the origin of the 
mapping referenced to the airfoil 
leading edge. Recommend approximately 
X(LE) +1/2 leading edge radius. 


• 

YSING. - 

Y coordinate of the origin of the 
mapping referenced to the airfoil 
leading edge. Recommend approximately 
Y(LE). 

16 

1 

DESC. 

Description for cards in Read Order 17 
Format ( 8 A 10 ). 


17 MU XP(I), YP(I) 

Format ( 8 EIO. 7 ) 

XP(l). - X coordinate of airfoil upper surface, 
ordered leading edge to trailing edge, 

YP(l). - Y coordinate of airfoil upper surface, 
ordered leading edge to trailing edge. 
Note that there is only one pair of 
coordinates per card. 

If airfoil section is not symmetric (FSYM<1,) the airfoil lower surface 
coordinates must be read here. For symmetric airfoil (FSYM^l,), skip the 
two Read Orders I 8 and 19* 
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Read Number 

Order Cards Description and Comments 


18 

1 

DESC . ~ 

Description for 
Format (8A10). 

cards in Read Order 19 

19 

FNL 

VAL, DUM 

- Format (8E10.7) 





VAL. - 

X coordinate of 
ordered leading 

airfoil 
edge to 

lower surface, 
trailing edge. 



DUM. - 

Y coordinate of 
ordered leading 
Note that there 
coordinates per 

airfoil 
edge to 
is only 
card. 

lower surface, 
trailing edge, 
one pair of 


Read Orders 10 through 19 complete the input for one span station. As 
indicated above Read Order 8, at least Read Orders 10 and 11 must be repeated 
for the remaining FNC-1 sections vhen FNC > 2. 


The above Read Orders complete the input for a computational case or a 
sequence of cases if obtained by having FNC < 2. Additional cases requiring 

different input are obtained by repeating the above Read Orders. Recommend that 
cases be run one at a time. 


The program terminates by reading the first three Read Orders with FNX<l,j 
that is, last three cards for a normal stop should be: 


1 

1 

TITLE . 

- End of Calculation 

2 

1 

DESC. 

- Description for card in Read Order 3 

3 

1 


0* • . • • 


SAMPLE CASE 


A simple wing, shown in Figure 1, is used here for the sample case; more 
realistic wing shapes are considered in references 1 and 2. This wing has 
a rect^gular planform of aspect ratio 32/9 with a NACA 65AOO6 airfoil section 
(ref. 8) at all span stations- The free stream, at Mach number .9> is at 2^ 
incidence and 0^ yaw with respect to the wing. Note that the freestream 
direction is consistantly indicated by a large open arrow on the figures. 

Figure 2 depicts how the wing section planes are transformed from 
physical space to computational space for the present sample case. In 
Figirre 3» a portion of the equally spaced computational grid in the wing 
section plane is plotted in physical space. This plot, however, was obtained 


8 



from a separate program,; supplied "by J. D. Keller of NASA Langley Research 
Center. . ^ ^ 

Appendix B is a listing of the; input cards for the sample case;. Read 
Order notations on the right side identify these cards with the input 
description of the last section. 

Appendix C is a collection of figures which show some typical output 
results for the sample case. Note that FPLOT was 2 in the input (Appendix B) 
so there are two types of plots; those from THREED and those from GRAPH. 

The bottom line of printing is common to all plots and contain^ the 
NASA LRC computer process number, date and time of run. Another line common 
to all plots contains the freestream Mach number (m), yaw single (YAW), 
and incidence angle (ALF). Figures Cl, C2, and C3 show the plots obtained 
from subroutine THREED. The first shows interpolated airfoil sections at 
each of the spanwise computational planes for the fine grid. (For this 
sample case, all sections are identical.) This plot serves to locate the 
CP origins on the second and third plots, where curves for all sections 
are superimposed but shifted with respect to one another. Figure C2 shows 
the chordwise distributions of (negative) surface pressure coefficient on 
the upper surface of the wing at all spanwise coi^putational planes of the fine 
grid, whereas Figure C3 gives that for the wing lower surface. The sharp 
downward breaks in the cxirves (at about 30-60% chord) are shock waves. On 
these plots L/D, CL and CD are for the entire wing; CL is the inviscid 
lift coefficient, CD is the drag coefficient (composed of that obtained from 
integrating the inviscid pressures plus that read in as input (CDO) due to 
skin friction, and L/D is the ratio CL/CD. A set of plots like these is 
obtained on each grid refinement; for the sample case there were two grids so 
two sets were obtained. 

Figures Ck and C5 are typical of the plots obtained from subroutine 
GRAPH. In these, the chordwise distribution of surface pressure coefficient 
for both the upper and lower wing surfaces is shown on one figure, a 
separate figure for each spanwise conputational plane. The value of Z, 
shown on the next to last line at the bottom of the plot, is the spdnwise 
coordinate in the physical plane (See Fig. l) and thus identifies the 
section location. The values of CL and CD shown on these plots are the 
inviscid section lift and drag coefficients, respectively. Figure Gh 
is for the root section (Z = 0.0) ^ile Figure C5 is for the tip station 
(Z = l60.0). The long tick mark on the ordinate is the sonic value of 
pressure coefficient. A set of plots (one for each spanwise station) is 
obtained for each grid refinement. * 

Figures C6, CT, and C8 sure charts which show the local Mach mnnljer 
(multiplied by 100) in various computational planes. Figxire C6 is at the 
root section of the wing (Z = 0.0 ); minus signs on it denote the lender half 
plane. A solid line outlines that portion of computational grid shown in 
Figure 3. The supersonic bubbles on both upper and lower wing sxirfaces are 
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clearly visible. Figures C7 and C8 show the Mach charts in the wing plane 
at the upper and lower STirfaces, respectively. The wing planform (in the 
con^jutational plane) is outlined on each figure and again the supersonic 
bubbles on both surfaces of the wing are clearly visible. One set of these 
charts is made on each grid refinement. Those for the fine grid of the sample 
case have been given in figures C7 and C8. 
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0. NACA 63AOO6 90. 

section 

Figure 1 .- Geometry of ving for sample case. 
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APPEISTDIX A - Langley Computer Program A4231 


NASA- LANGLEY RESEARCH CENTER 


01 i 

LAR 


7 rttoi'.KAM r«0 


A4231 


COMPUTER PROGRAM ABSTRACT 


01 14 0AT€ 

081073 


01 ro title of phowau fsi cmaracteks maximum^ 


Transonic Flov Analysis-3D airfoils/Jameson (parabolic) 



07 77 

language 

NO. 1 

or :I7 

language 

NO. 7 

m 

F^R6H 


WHO 

os 14 CONTACT 


PAYiEin PROGRAM 


07 14 

CATEGORY 


07 IS SITE 

NYU 


07 IB PROGRAM NO 


07 V KEY words I BMAXMUM separated BY COMMAS! 

transonic flow, embedded shocks, Isentroplc supersonic flow. 


subsonic flow 


WHO TO CONTACT ABOUT THE PROGRAM 


RMDavis 


ns 74 SITE 


LAR 


Wr, T1 or<;n code 


33.300 


TiS J7 PROItCT NO 


R2564 


tTj 4S 

NASA 

CCNUR 


4« status 


CD UNDER nCVf.LOm.NT 
0. operational 
C3 c. COMPLETED 


ns *3 

r~1 * This ('RfKiRAM 


IS NO r t OR 
SHARING 


OATtS 


05 sa WCVIStONCOOt 


TUtt AND COST FOR OtVEt OPMFHT 


OS wihitiated 


0873 


as S4 cOMriETp.D 


□ A RtViSlON 
O 0 CANCELLATION 


OS SO MANMONTHS 


>9 60 61 62 63 


OS 64 MACHINE 
HOURS 


64 6S 66 67 66| 


OS '69 COMPUTER TYPE 


6000 


OS 14 TOTAL COST 
iDOLLARSi 


k k .1 I A 1 


74 ;S ?6 17 74 79 M 


ELITE MARGIN 


CARO NUMBER 


ABSTRACT 


This program performs calculations for Inviscid isen^roDlc 


transonic flow over 3D airfoils with straight ('may be ya^ 


1 )11^ not gvepi^) leading edges. The fre 6 ’^stream Mach numbBr Is 


rttfiH only the laentroplc aasumptlciia.. Weak , 


arft iQcftted vhere ever they Qccur In, th^_0_ow^ 


The f 1 n1 terrdl ff finance form nf thfi fiiTT fiquatlxai for the TelQclt^ 


pot e ntial i e eo lv e d by the m e th od o f rel aYe t t lnn, after the 


flow e xt e ylop to th e airf o il le m a pp e d t o th e iipp e r h a lf 


piano i — Tho mapping — ppooodure allow e e a c a e t ga tl e f a at i cai o f . th e} 


boundary e ondltiono and uoo of oup e poonlo fr ee s tg e am v e looltl ec « 

tm — A.- Ajn L. 


finit e diff ere n ce oper ato r la ^lo c ally rotated* — 4» 


a u p e r a e nt e flo w reg ion s sa as to properly aco e uat for th e 


d o main o f d epe nd e n ce . — Th e re la a eation al ger l t hB has b e en 


sta b ilized usin g c rite r ia f ro m a tiawlik e anal o g y . 


Figure Al.- Coniputer program abstract. 
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APPENDIX A - Langley Computer Program 


GEOM 

COORD 

SURF 

ESTIM 

REFIN 

INTPL 

SPLIF 

MKFLO 

YSWEEP 

XSWEEP 

VELO 

FORCF 

CPLOT 

TOTFOR 

CHARTZ 

CHARTY 

GRAPH 

THREED 

PLOT * 

SYMBOL * 

AXIS * 

LINE * 


reads geometric definition of wing 

sets up stretched parabolic and spanwise coordinates 

interpolates mapped wing surface at computational mesh points 

computes initial estimate of reduced potential 

halves computational mesh size 

interpolates using Taylor series 

performs cubic spline fit 

solves mixed- flow equations (rotated difference scheme) 
relaxes horizontal lines, sweeping in y-direction 
relaxes vertical lines, sweeping in x-direction 
calculates surface velocity 
calculates section force coefficients 

generates printer plots of CP at equal computational intervals 

calculates total force coefficients 

generates Mach number chajrt in plane of wing section 

generates Mach number chart in plane of wing planform 

generates Varlan plots 

generates three-dimensional plots 

moves pen (call CALPLT) 

plots symbols (call NOTATE) 

draws axes (call AXES) 

plots arrays 


■K expedient modifications were made so that NYU plotting subroutines were 
compatible with NASA LRC CDC system. 


Figure A5. - Function of subroutines 



APPENDIX B - INPUT FOR SAMPLE CASE 








Read 

Order 

RFCTANGULAR W I NG - 

NACA 63A006 AR-32/9 



1 

FNX 

FNY 

FN2 

FPLOT 



2 

48. 

B « 

16* 

2. . „ . 




FIT' 

COVO 

PI 0 

P20 P30 

0ETAO STR1P0 

FMALF 

U 

' zzo\ 



1 .6 

1 • I * 

07 1 . 

t • 

5 

1 1 0. “ 

'.O'OWl 

1 .6 

1 . 1 . 

0. 1 • 

0. 


■ FMACH 


AL 

COO 



6 

■ * 9 

0.“ ^ 

2. 

0. 



7 

FNC 






8 

?. 






9 


CHORD 

TNICkS 

“AL FSEC 



10 ^ 

11 ! 

-160. 

’90. 

1 . 

0* 1 * 



F5YM 

FNU 

FNL 




12 ! 

1 • 

36* 

36. 




13 1 

TRL 

SLT 

XSING 

YStNG 



Ih I 

7.0770 

' 0. 

-.1325 ■■ 

0. 



15 

XP’ 

VP 


1 UPPER SURFACE 1 



i6 : 

0, 

0. 





I 

.05555? 

.16235 





j 

• Ill 11 

,22444 





• 

.2222? 

.32222 





] 

..13333 

.401 1 1 i 






'•4444 4 

.46777“ 






.=^5555 

V52 " 






.63333 

.62277 






1 . Ti n I 

— ^.mii ■ 






1 .■666'66 

' .'364 70 






2.2222 

' '".98961 






"2Vr7T7 







■ 3.3333 

1 « 1 962 





t 

1 

3 .8688 

1 .2877 
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APPENDIX C - Some Output Results For San^le Case 



Note : Spanwise variable 
is not scaled to 
chordwise variable in 
this plot* Instead, the 
wing sections shown 
here located the Cp 
origins in the following 
plots* 


VIEW or WING 

RECTANGULAR WING - NfiCA 63R006 SECTION AR=32/3 
M = .900 YAW = 0.00 RLF = 2.00 

L/D = 23.05 CL = .1876 CD = .0081 

LB71564 09/13/73 07.37.22. 


Figure Cl.- View of wing showing interpolated edrfoil sections on fine grid. 
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APPENDIX C - Some Output Results For Sample Case 



UPPER SURFACE PRESSURE 

RECTANGULAR WING - NBCA 63A006 SECTION AR=32/3 
M = -300 YAW = 0.00 ALE = 2-00 

U/D = 23.05 CL = .1876 CD = .0081 

LB71564 09/13/73 07.37.22. 


Figure C2*- Chordwise distributions of upper-- surface wing pressure 

coefficients at spanwise computational planes of fine grid, 
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APPENDIX C - Some Output Results For Sample Case 



LOUER SURFRCE PRESSURE 

RECTRNGULfIR WCN6 - NflCfl 63R006 SECTION flR=32/3 
M = .300 YflW = 0.00 RLF - 2.00 

L/D = 23.05 CL = .1876 CD = .0081 

LB7156H 03/13/73 07.37.22. 


Figure C3.- Chordwise distributions of lover- surface ving pressure 

coefficients at spanwise computational planes of fine grid* 
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■ .. APPfflDIX C - Some Output Results For Sample Case 


+ upper wing surface 
X lover ving surface 



RECTRNGULflR WING - NRCR 63R006 SECTION RR-32/3 
M .= .300 YRW = 0.00 RLF = 2.00 

Z= 0-00 CL = .2261 CD = .0107 

LB71564 03/13/73 07.37.22. 


Figure Chordwise distributions of pressure coefficient at root 

section of wing* 
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APPENDIX G - Some Output Results For Sample Case 



RECTflNGULFIR WING - NRCfl 63R006 SECTION BR=32/3 
M = .300 YRW = 0.00 RLF = 2.00 

Z = 160.00 CL = .0306 CD = .0036 

LB71564 03/13/73 07.37.22. 

Figure Chordwise distribution of pressure coefficients at 

tip section of wing* 
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APPEHDIX C - Some Output Results For Sample Case 
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Figure C6.- Mach chart at root section of wing in computation plane. 
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APPENDIX C - Some Output Results For Sample Case 


Free 


stream 



SUBfiACP MACH NO CHART TN «1NC PLANf 



PQ 

H9 

89 

89 

09 

05 

B9 

P*= 

85 

6® 

09 

B9 

89 

0® 

89 

09 

89 

09 

89 

09 

89 

09 

89 

89 

09 

89 

89 

09 

09 

89 

09 




89 

89 

89 

89 

39 

0* 

B9 

99 

85 

89 

09 

89 

89 

09 

«V 

89 

b9 

89 

89 

59 

09 

B9 

89 

09 

09 

R9 

09 

89 

89 

09 



89 

89 

B9 

BP 

85 

0P 

BP 

H9 

89 

89 

89 

05 

fi® 

89 

89 

89 

89 

85 

89 

89 

6® 

09 

HO 

B® 

09 

B9 

89 

89 

09 

89 


fj«J 

89 

89 

B9 

69 

39 

89 

89 

88 

68 

88 

BB 

39 

83 

88 

8b 

88 

bB 

80 

08 

00 

00 

80 

69 

89 

89 

69 

09 

09 

89 

89 



89 

89 

89 

66 

08 

88 

88 

8B 

00 

87 

67 

07 

87 

07 

87 

87 

07 

87 

07 

07 

86 

00 

80 

80 

88 

80 

H9 

89 

09 

09 


flS 

89 

89 

80 

80 

08 

87 

87 

07 

B6 

06 

06 

86 

06 

86 

8b 

86 

06 

86 

86 

86 

66 

87 

07 

07 

80 

B0 

00 

89 

B9 

89 


69 

09 

89 

63 

0T 

87 

06 

86 

95 

B* 

35 

85 

05 

85 

04 

84 

»4 

04 

B5 

85 

85 

05 

05 

66 

P6 

87 

87 

88 

89 

09 

89 



90 

09 

89 

88 

87 

06 

85 

B4 

64 

8? 

B3 

e? 

03 

83 

63 

85 

03 

03 

63 

Bl 

03 

83 

b4 

64 

05 

86 

07 

88 

89 

S9 

90 


90 

89 

09 

88 

85 

04 

83 

0? 

H? 

81 

61 

81 

0! 

81 

Si 

81 

01 

0» 

Bt 

81 

Bl 

61 

B2 

82 

6? 

P4 

85 

00 

89 

H9 

90 


90 

89 

89 

86 

B6 

•02 

8} 

00 

75 

75 

79 

79 

79 

79 

7b 

78 

7b 

79 

79 

79 

79 

79 

79 

60 

01 

83 

86 

00 

B9 

89 

90 


90 

90 

B9 

06 

85 

60 

77 

77 

76 

76 

76 

76 

7e 

76 

76 

Tb 

76 

76 

76 

76 

76 

76 

76 

77 

77 

80 

05 

88 

89 

90 

00 


90 

90 

69 

88 

85 

76 

73 

7? 

72 

72 

72 

72 

72 

72 

7i 

72 

72 

72 

72 

72 

72 

72 

72 

73 

73 

76 

05 

80 

89 

90 

90 


90 

90 

89 

89 

66 

70 

60 

67 

67 

67 

67 

67 

67 

6*^ 

67 

67 

67 

67 

67 

67 

67 

47 

67 

67 

63 

70 

86 

09 

89 

90 

90 



90 

90 

OO 

8® 

86 

62 

6C 

tC 

59 

59 

59 

59 

59 

59 

59 

59 

. 55 

. 5^ 

. 59 

59 

59 

59 

59 

60 

60 

62 

86 

09 

91 

90 

90 


90 

90 

90 

89 


47 

46 

46 

46 

46 

46 

46 

46 

46, 

^9 

46 

46 

46 

4£ 

46 

46 

46 

46 

46 

46 

47 

87 

89 

90 

90 

90 


90 

90 

90 

89 

0T ] 

19 

21 

72 

•j 3 

?? 

24 

24 

24 

24 

24 

24 

24 

24 

74 

24 

>4 

23 

?2 

22 

21 

19 

87 

89 

90 

90 

90 


90 

90 

90 

89 

07 t 

78 

80 

81 

B2 

83 

03 

03 

B4 

84 

84 

84 

04 

04 

84 

83 

03 

03 

07 

Bl 

BO 

78 

BT 

89 

90 

90 

90 


90 

90 

90 

09 

87 IOC 

101 

1C2 


10? 

103 

104 

1^4 

104 

104 

104 

104 

104 

104 

104 

1 0? 

1?3 

103 

102 

19’ 

100 

87 

09 

90 

90 

OQ 


90 

90 

90 

69 

86 104 

1C3 

1C4 

1C4 

105 

! 35 

105 

105 

106 

106 

106 

106 

106 

105 

105 

105 

105 

104 

104 

102 

104 

68 

69 

90 

90 

90 


90 

90 

90 

89 

80 »05 

104 

1 C4 

? 04 

104 

105 

!05 

105 

105 

105 

105 

105 

105 

105 

10 5 

’ 05 

104 

104 

104 

1 04 

105 

00 

89 

91 

6Q 

60 


90 

90 

90 

90 

89 106 

lot 

!C6 

IC6 

!06 

106 

106 

1 C6 

1 06 

106 

106 

106 

106 

1C6 

106 

106 

106 

106 

106 

106 

lot 

89 

90 

91 

90 

90 


90 

90 

90 

9^ 

99 107 

107 

1C7 

»C7 

107 

107 

107 

’ 07 

lOT 

107 

1J7 

1W7 

107 

107 

107 

107 

107 

107 

107 

107 

107 

90 

90 

90 

90 

9-0 


90 

OO 

oo 

90 

91 1C7 

108 

1 05 

) 03 

’06 

100 

108 

100 

108 

108 

108 

lUe 

1U0 

!0B 

!0B 

lOS 

108 

1 Ofl 

108 

108 

107 

91 

90 

90 

90 

90 



90 

90 

91 

92 ’ 

!07 

109 

?09 

7 C9 

109 

109 

109 

109 

109 

109 

109 

109 

1 09 

1C9 

l'>9 

109 

1 '■•9 

109 

109 

109 

107 

9? 

91 

90 

90 

00 


90 

OQ 

91 

91 

9? 

< 10 

ll'* 

!!*“ 

no 

no 

no 

no 

in 

111 

111 

111 

111 

no 

! 10 


no 

no 

no 

IID 

\fi6 

93 

91 

9! 

90 

90 



50 

90 

91 

9? 

93 105 

110 

Ul 

112 

112 

?12 

112 

11? 

L17 

112 

112 

112 

112 

*12 

112 

112 

112 

112 

in 

no 

t 

92 

92 

91 

90 

90 


50 

91 

91 

9? 

94 105 

lie 

1 12 

m 

113 

! 7 7 

U? 

L! 

113 

113 

113 

113 

112 

1 12 

LL3 

113 

113 

in 

112 

no 

105 

94 

92 

91 

91 

90 


90 

91 

91 

93 

95 104 

110 

ill 

U4 

114 

114 

114 

ii4 

114 

114 

114 

114 

114 


il4 

114 

1 1 4 . 

114 

U3 

no 

105 

95 

93 

91 

9! 

90 


90 

9\ 


9? 

96 104 

110 

n? 

7 ! 5 

115 

!15 

115 

115 

115 

115 

115 

U5 

115 

115 

1? 5 

J15 

1 15 

1*5 

m 

no 

105 

96 

93 

92 

91 

90 


90 

9* 

9? 

94 

96 105 

110 

1 14 

IIS 

U6 

116 

116 

116 

116 

lib 

lib 

Ub 

U6 

U6 

116 

116 

116 

115 

U4 

110 

’05 

96 

94 

93 

91 

90 


90 

91 

9? 

94 

97 105 

110 

714 

116 

117 

7.18 

no 

170 

lie 

118 

118 

lib 

110 

lie 

ua 

118 

T17 

116 

114 

I’O 

105 

97 

94 

93 

9’ 

90 


_9C 

91 

93 

95 

9P 106 

111 

U5 

l!7 

119 

»19 

119 

7 15 

1 19 


.119 

119. 

119 

U9_ 

.119 

119 

119 

U7 

ns 

m 

106 

9B 

95 

93 

91 

90 


90 

91 

9? 

95 

98 106 

m 

115 

118 

1 20 

121 

121 

7.21 

121 

121 

121 

121 

121 

121 

1 21 

12! 

120 

118 

115 

112 

106 

69 

95 

93 

91 

90 


90 

91 

93 

95 

99 107 

112 

116 

! 79 

12J 

122 

123 

123 

123 

123 

123 

123 

123 

l?3 

123 

122 

121 

119 

116 

112 

107 

99 

95 

93 

91 

90 



90 

91 

93 

95 

99 if>7 

113 

116 

120 


134 

124 

124 

125 

125 

12 s 

125 

125 

?24 

124 

124 

122 

12'’ 

117 

11 ? 

1?7 

99 

95 

01 

®1 

0(5 


90 

9» 

93 

. 95.- 

98 !07 

1 12 

1»7 

! 20 

123 

1.25 

126 

126 

126 

126 

12b 

12b 

126 

126 

126 

125 

123 

120 

117 

112 

10 7 

93 

95 

93 

c\ 

90 


90 

91 

93 

94 

97 ! 

,05 

lie 

1!4 

119 

7.22 

125 

’27 

127 

1?B 

128 

128 

120 

i28 

127 

1 27 

125 

122 

119 

U5 

no 

105. 

6T 

05 

93 

91 

00 



90 

°.l 



96 ! 

101 

103 

JC6 

1C9 

113. 

..119 

12? 

125 

126 

127 

127 

127 

126 

125. 

123 

119 

114 

109 

:o6 

-102. 

lOl 

96 

94 

93 

91 

90 


_ 90 _ 

91 

53 

93 

95. 

90 

.97. 

. 

58 

100 

107 

105 

107 

109 

no 

111 

lie 

IC9 

!07 

105 

102 

100 

98 

97 

97 

98 

95 

93 

92 

91 

90 



-^^-1 

91 

91 ■ 


93 

95 

45 

94 

94 

93 

92 

9? 


'94 

94 

95 

94 

54 

93 

92 

92 

93 

94 

94 

95 

95 

93 

93 

91 

91 

90 


50 


91 

91 

53 

93 

92 



91 

.91 

90 

9».. 

97 

91 

91 

91 

_5l. 

-91. 

90 

*>'. 

. 91 

91 

91 

9?_ 

93 

92 

91 

91 

91 

90 




_90. 

C5 

90 

91 

91 

90 

0P 

89. 

..98 

86 

88 

88 

66 

88 

- SB 

88 

8B 

86 

88 

08. 

_.86. 

. 89 

.09 

.90 

91 

91 

90 

90_ 

90 

90 



90 

~90' 

50 

90 

90 

89 

07 

86 

Be' 

86 

86 

B6 

66 

06 

86 

bb 

8b 

06 

86 

66 

66 

e6 

06 

86 

87, 

89 

90 

90 

90 

90 

90 




90 

_?0 8? 

BB 

65 

,-.07. 

__ a? 

a? 

82 

82 

82 

82 

82 


02 

_S2. 


82 

82. 

82 

82 

02 

8? 

83 

85 

93 

S9 

90 

90 

90 



90 " 

90 

.89 

09 

8« 


81 

"so 

eo 

80 

ao 

00 

' 8 0 ' 

'bo 

3C 

ap 

80 

80 

80 

80 

80 

ao 

80 

00 

8t 

82 

. 8 8 

89 

B9 

90 

®0 



1 9C_I 

89 B9 

__?A. 

.00 

84 

82 

82 

82 

82 

81 

31 

0J 

8! 

ai 

81 

bl 

. -01 . 


01 

.81 

. 62_ 

_S2 

..92 

,_e? 

. 04 

..88 

89 

39 

90 

OQ 


90 

89 

89 

89 

89 

86 

B5 

85 

B5 

84 

94 

04 

64 

84 

84 

84 

84 

84 

B4 

04 

84 

84 

85 

85 

85 

B6 

80 

89 

89 

69 

90 


90 

89 

89 

89 

66 

87 

87 

86 

86 

06 

86 

85 

85 

85 

85 

«5 

85 

05 

85 

05 

66 

86 

86 

36 

B7 

87 

88 

89 

89 

89 

®0 


90 

89 

89 

69 

89 

68 

87 

87 

0T 

87 

86 

06 

86 

86 

8o 

86 

86 

66 

86 

06 

65 

07 

07 

87 

07 

88 

89 

89 

89 

89 

90 


89 

89 

89 

89 

69 

66 

08 

89 

87 

87 

97 

07 

07 

97 

87 

87 

87 

87 

87 

07 

87 

87 

87 

88 

86 

88 

89 

89 

89 

89 

09 


0O 

89 

89 

09 

84 

89 

88 

60 

00 

38 

B8 

08 

88 

87 

07 

87 

87 

07 

88 

88 

B3 

B8 

88 

86 

BS 

89 

89 

09 

89 

09 

89 


69 

89 

69 

B9 

69 

B9 

89 

8B 

68 

88 

89 

BB 

88 

85 

08 

88 

8b 

88 

66 

88 

88 

B8 

88 

88 

39 

89 

39 

89 

89 

89 

89 


89. _ 

89 

89 

09 . 

.09 

85 

85 

. 89 

08 

00 

. 80 

08 

88 

aa 

88 

68 

. 88 

. 80. 

- 88. 

88 

63 

B0 

80 

89 

89 

. B9 

89 

89 

89 

89 

89 


99 

89 

B9 

B9 

90 

89 

89 

89 

65 

85 

89 

09 

89 

SB 

68 

88 

8b 

88 

89 

89 

89 

99 

89 

39 

89 

89 

90 

69 

89 

89 

69 


89 

89 

89 

8P 

9 0 

39 

89 

fic 

89 

89 

89 

09 

09 

89 

89 

89 

89 

B9 

89 

89 

89 

89 

89 

09 

09 

89 

90 

89 

09 

89 

09 


09 

09 

09 

89 

90 

09 

89 

85 

09 

89 

89 

89 

89 

B9 

89 

89 . 

89 

89 

89 

89 

89 

.09 

89 

89 

89 

89 

90 

.89 

89 

89 

09 


89 

89 

P9 

89 

90 

09 

89 

89 

09 

89 

89 

09 

89 

99 

89 

B9 

B9 

09 

89 

89 

89 

09 

89 

. 89 

89 

89 

90 

39 

09 

89 

89 


89 

89 

09 

.90 

93 

90 

89 

89 

89 

89 

89 

09 

89 

89 

89 

89 

89 

89 

89 

89 

89 

89 

89 

39 

89 

90 

90 

90 

89 

09 

89 

-- 

89 

89 

09 

90 

9C 

90 

89 

8® 

8® 

99 

69 

99 

89 

B« 

89 

«9. 

,05 

89 

_ gs 

. o9 

89 

89 

8® 

._Sv 

-89. 

- 90. 

.. 90 

90 

09 

09 

89 


89 

89 

90 

90 

9C 

90 

89 

89 

89 

89 

89 

B9 

89 

89 

89 

89 

09 

09 

89 

89 

89 

89 

89 

89 

89 

90 

90 

90 

90 

89 

89 


69 

90 

50 

90 

9C 

90 

90 

89 

89 

89 

89 

89 

89 

69 

89 

85 

89 

89 

09 

89 

09 

89 

69 

09 

89 

90 

90 

90 

90 

90 

89 


90 

90 

90 

90 

9C 

9C 

90 

89 

89 

89 

99 

89 

89 

09 

89 

89 

09 

09 

89 

89 

89 

89 

89 

89 

90 

90 

90 

90 

99 

90 

90 


90 

90 

50 

90 

90 

89 

89 

89. 

09 

39 

99 

B9 

89 

89 

09 

89 

-09 

.. 09 

89 

69 

89 

80 

69 

89 

89 

. 89_ 

90 

90 

90 

90 

90 








— ■ 


— 







— 

--- - 

— 

— 




— 








, 


_ 






Figure CT*^ Mach chart at upper surface of ving in computational plane. 


2k 



APPENDIX C - Some Output Results For Sample Case 
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Figure C8.- Mach chart at lower surface of wing in 
computational plane. 
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